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Abstract 



Prom a general model of the Mn oxides Ri- x A x MnO^, we derive an effective 
Hamiltonian in the low-energy subspace using the projection operator method, 
in which a novel coupling between the spin and orbital degrees of freedom is 
included. A phase diagram for temperature T versus doping concentration x is 
computed by means of Monte Carlo simulation. Our result is consistent with 
experimental observations in the Mn oxides with relatively wide conduction 
band, such as Pri_ x Sr x Mn03 and Lai_ x Sr x Mn03. According to the obtained 
orbital ordering, we also predict that the motion of charge carriers transforms 
from three-dimensional to two-dimensional as x is increased beyond a critical 
value. 
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The peculiar magnetic and transport properties including the colossal magnetoresistance 
(CMR) of the manganites Ri_ x A x Mn0 3 (where R =La, Nd or Pr and A =Ca, Sr, Ba or Pb) 
have become the main focus of many recent research activities. In undoped -RMnOs, each 
Mn atom has four outshell 3c? electrons, three localized t2 9 electrons forming an S = 3/2 core 
spin, and one itinerant electron filling into two degenerate e g orbitals, namely, |+) = d 3z 2_ T 2 
and |— ) = d x 2_ y 2. i?Mn03 is usually an A-type antiferromagnetic (A-AF) insulator [|l],@, 
where the local spins are aligned ferromagnetically in layers parallel to the ab plane and 
antiferromagnetically along the c axis. Hole doping by partial substitution of R 3+ atoms by 
A 2+ atoms soon destroys the A-AF order, and strong ferromagnetism occurs for 0.1 ~ x ~ 

Near x = 0.5, the manganites exhibit different magnetic structures, as controlled 
by the conduction band width W. some manganites with relatively large W, such as 
Pr .5Sr .5MnO 3 |||7||i , are in an A-AF state below a critical temperature T N and changes 
to a ferromagnetic (FM) state above through a discontinuous phase transition. Some 
manganites with small W, e.g., Pro.sCao.sMnOs || and Lao.5Cao.5Mn03 |10" |, exhibit a more 



complicated CE-type antiferromagnetic (CE-AF) order with charge ordering at low tem- 
peratures, where the spins arrange ferromagnetically in zig-zag chains and antiferromagnet- 
ically between neighboring chains. In some other manganites with intermediate W, such 
as Ndi-^Sr^MnOa |j8,nfl, CE-AF and A-AF orders may coexist in the systems, and their 



volume fractions change rapidly with small change in x. 

For further doping, the A-AF state appears to be the fundamental ground state spin 
configuration in many manganites. It has been observed in Lai-^Sr^MnOa with 0.52 < x < 
0.58 0, Ndi-a-Sr^MnOa with 0.53 < x < 0.62 |TJ and Pn^Sr^MnOg with 0.48 < x < 



0.6 II . Near the ending doping concentration x = 1.0, the manganites develop an isotropic 
G-type antiferromagnetic (G-AF) ground state []IJ. In the region between x ~ 0.6 and the 
ending concentration, most manganites are non- ferromagnetic insulators [§,|l2j) but their 
microscopic spin structures have so far been seldom studied in experiments. 



The conventional double exchange (DE) model |T3Hl5|] for the Mn oxides usually pre- 
dicts strong ferromagnetism symmetrically about x = 0.5, which does not account for the 
antiferromagnetic phase at x ~ 0.5 observed experimentally. While Jahn- Teller (J-T) ef- 
fect [[I6| , |r7|1 , orbital degrees of freedom |Tj|, and Coulomb interactions [|5"1,|18| have been 



proposed as necessary extensions of the DE model, theoretical phase diagram in the T ~ x 
plane commensurating with experimental observations has yet to be obtained. 

In this work, starting from a general model for the manganites we derive a modified 
DE Hamiltonian, which by comparison with the conventional DE model contains an orbital 
dependence in the hopping integral and a novel effective coupling between the spin and 
orbital degrees of freedom originating from the strong correlations. The magnetic and orbital 
phase diagram of this Hamiltonian is calculated using Monte Carlo simulation and simulated 
annealing technique [[HJ. Our results obtained are in good agreement with the experimental 



measurements in the manganites with relatively wide conduction band in the entire T ~ x 
plane. 

With the two-fold degeneracy of the e g orbitals included, the total Hamiltonian for the 
Mn oxides can be written as 

H = — a \as a ji3s — K cr ss i ■ Sia\ as ai as i + U ^ n^n^ 

ijaj3 iass' ia 
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j _ k 

+ U'J2 n i+ sni-s> +y£scs r s^ T a/3 ■ Qia\ as a i0s + - ^ Q 2 . (1) 

iss' ij iaf3s i 

Here, aj as annihilates an electron at site i in e g orbital a(= +, or — ) with spin s(=|, or 
|). K is the Hund's rule coupling between e g electrons and localized spins, where er are 
the Pauli matrices in the e g spin space and Sj = Sj/S 1 a unit vector in the direction of the 
classical localized spin Sj. U and U' represent the intraorbital and interorbital Coulomb 
repulsions, respectively. J s stands for the superexchange coupling of the localized t 2g spins. 
g is the J-T coupling between the e g electrons and local J-T lattice distortions Qj, where r 
are the Pauli matrices in the orbital space and Qj = (Qi X ,Qiz) describes the two e g modes 
of the J-T distortions. The orbital degrees of freedom are also called the isospin of the e g 
electrons, for they are similar to the spin degrees of freedom. The orbital-dependent hopping 
integrals |L8| tff are elements of the 2x2 matrix = t(l + r ■ ny). Here, % = 11,5 with 



5 = x,y, and z for hopping along the x, y and z direction, where n^, = (— v / 3/2,0, —1/2), 
n sz = (a/3/2,0, —1/2) and n z = (0,0, 1) are three unit vectors distributed symmetrically in 
the x — z plane. Introducing the three unit vectors is a key step in our theory, which will 
allow us to express the effective Hamiltonian in a compact and transparent form. 

In the case of large U, U', K and J-T coupling g, we can derive an effective Hamiltonian 
in the low-energy subspace by following a similar procedure to the derivation of the t — J 



model from the Hubbard model p0[ . The total Fock space is divided into two subspaces: 
the low-energy subspace where each Mn site is occupied by at most one e g electron and 
the spin and isospin of the e g electron are parallel to local s$ and Qj, respectively, and 
the reminder subspace, namely, the high-energy subspace. Using the projection operator 
method and treating the kinetic energy term as a perturbation P0,|2"T|, we obtain an effective 
Hamiltonian in the low-energy subspace 

H eff = ~ tij4 d j + y XX 1 + W»i)si • s i 
J' 

- 17 X^ n i(l + As * ■ s i)[l - ( r * ■ n ij)( T j ■ n ij)} . ( 2 ) 

Z ij 

where rii = djdi, and di annihilates an e g electron with spin and isospin parallel to Sj and Qj. 
Here, 7 = 2/ 3 /J s , J' = h + I 2 - h, and A = (h - h + h)/J\ where h = t 2 /(U' + 2n e Ej), 
I 2 = t 2 /(U' + 2n e Ej + IK) and J 3 = t 2 /(U + 2K) with Ej = g 2 /k as the J — T energy 
and Tig — 1 OC cLS the average electron density. The isospin and lattice distortion are forced 
mutually parallel by the strong J-T coupling. Tj is a unit vector describing the orientation 
of the isospin or which orbital state on Mn site % is occupied by an e g electron. If Tj makes 
an angle <pi with the z axis, then the e g electron at site i occupies the orbital state 

= cos(4>i/2)d 3z 2_ r 2 + sm((p i /2)d x 2_ y 2 , (3) 

and it also represents the direction of the lattice distortion; i.e, Tj = Qj/Qj. The effective 
hopping integral in Eq. (2) can be written as 



tij = t J (1 + Si • Sj)/2J(l + n ■ n^) (1 + Tj ■ rijj) . (4) 
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The first square root together with t is the well-known hopping integral in the conventional 
DE model ||13|| , and the last term describes the dependence on the orbital alignments. 

Equation (2) is a modified DE model. The main merit of this model is that it contains 
only weak couplings between classical fields Sj and Tj. To show its relevance to the Mn 
oxides, we shall calculate the magnetic and orbital phase diagram. Focusing on the spin and 
isospin dynamics, we replace the charge operators in Eq. (2) by their averages rii ~ n e and 
d\dj ~ (d\dj). We further approximate (djdj) ~ (d\dj)o with (d\dj)o as the average of d\dj 
in a tight binding model. The resulting spin and isospin Hamiltonian is classic and can be 
studied by using Monte Carlo simulations together with simulated annealing technique [ITS 



We start from a sufficiently high temperature /c^T = 0.2t and work on a 14 x 14 x 14 cubic 
lattice. The temperature is decreased by 8% each step. For each step, about 2 x 10 6 spin and 
isospin configurations are sampled, and thermal average is calculated after initial sampling 
of 5 x 10 5 configurations. 

We consider first the undoped case, where n e = 1 or x = and the kinetic energy term 
in Eq. (2) vanishes. The ground state depends on only two dimensionless parameters A and 
J' I J where J — Js(l + 7). In the Monte Carlo simulation, as the temperature is decreased 
to a sufficiently low value, e.g., ksT = 10~ 2 J, the spins and isospins will eventually develop 
their ground state configurations, which can be examined directly. Numerical calculation 
indicates totally four possible ground states: G-AF, C-AF, A-AF and FM, where in the C-AF 
state the local spins arrange ferromagnetically in parallel chains and antiferromagnetically 
between neighboring chains. In all the four magnetic states, the e g electrons occupy two 
different orbital states on two sublattices, exhibiting alternating orbital order. For the FM 
and G-AF states, the e g electrons occupy an arbitrary pair of orthogonal orbital states in 
the two sublattices. For the A-AF state, the occupied orbital states in the two sublattices 
are (d 3z 2_ r 2 + d x i_ y -z) / \/2 and (d 3z 2_ r 2 — d x 2_ y 2)j \f 7 i. For the C-AF state, the two occupied 
orbital states are d 3x 2_ r 2 and d y 2_ z 2, if the ferromagnetic chains are assumed along the y 
direction. The obtained phase diagram in the A versus J' /J plane is given in Fig. 1. The 
phase transition between different phases with changing the parameters is of first order. 

In order to ensure an A-AF ground state for undoped i?Mn03, the parameters have to 
be chosen in the shaded region of Fig. 1. We will hereafter choose Js = 0.03t, U = 42t, 
U' = It, Ej = 3.5t and K = 9t as reasonable values for the Mn oxides, which correspond to 
the point indicated by the triangle in Fig. 1 for the undoped case. The phase diagram in the 
normalized temperature ksT/t versus doping concentration x plane is obtained numerically 
by calculating simultaneously various long-range order parameters such as the magnetization 
M and short-range correlation functions of the localized spins such as (sj • Sj) in three 
directions, as well as those of the isospins, as functions of ksT/t and x. The resulting phase 
diagram is shown in Fig. 2. 

In our phase diagram Fig. 2, there are an A-AF, a G-AF and a FM phase at low tem- 
peratures for x ~ 0, x ~ 1.0 and 0.1 ~ x ~ 0.4, respectively, which are consistent with 
experiments in most manganites . More importantly, we also observe an A-AF state for 
0.6 ~ x ~ 0.9 and a first-order FM to A-AF phase transition upon cooling around x ~ 0.5 
[as will also be seen in Fig. (4)]. This result agrees well with the experimental measurements 
in the manganites with relatively wide conduction band, e.g., Pri^Sr^MnOs It is 

also possible, by slightly tunning the input parameters, to make our phase diagram more 
comparable to that of Lai-zSr^MnOs, where a FM state exists at 0.1 ~ x ~ 0.5 B and an 
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A-AF state occurs for x > 0.52 112 



The magnetization M normalized by the saturation magnetization M s calculated using 
Monte Carlo technique is shown in Fig. 3 as a function of the normalized temperature 
for several values of doping concentration. For x = 0.1, 0.2 and 0.3, the magnetization 
approaches the saturation value M s at low temperatures, exhibiting strong ferromagnetism 
in this doping range. The magnetization for x = 0.4 drops sharply from a large value to zero 
at a low critical temperature, corresponding to a first-order FM to A-AF phase transition. 

Let us consider x = 0.5. In Fig. 4, the magnetization for different values of applied 
magnetic field is shown as a function of the normalized temperature. At zero magnetic field 
the magnetization first increases with decreasing temperature and then drops abruptly, after 
which the system enters an A-AF state. This magnetization curve resembles that measured 
in the Mn oxide Pro.sSro.sMnOa 0. With increasing the magnetic field, the magnetization 
increases obviously. For 2^bS*H /t = 0.024, which corresponds to H ~ 20T if t is taken 
to be 0.2eV as a reasonable value for the Mn oxides, the magnetization reaches about 30% 
even at very low temperatures. The sensitivity of the magnetization to the magnetic field 
may explain the sharp resistivity drop induced by applied magnetic field observed in the 
low-temperature A-AF phase of Pr .5Sr .5MnO3 0. 

Now we discuss the orbital structures in the phase diagram. As seen from the inset 
of Fig. 2, there are two different ground state orbital structures. In the region denoted 
by "antiparallel" , the e g electrons occupy two different orbital states (d 3z 2_ r 2 + d x 2_ y 2)/\/2 
and (dz z 2_ r 2 — d x 2_ y 2) / 'y2 in two sublattices, forming alternating orbital order. In the 
region denoted by "parallel", the e g electrons occupy the same orbital d x 2_ y 2 on each site, 
showing uniform orbital order. The spatial distribution of the charge density or the square 
of the orbital wave function in the two types of orbital orders is illustrated in Fig. 5. It is 
interesting to notice that the antiparallel orbital order can be reguarded approximately as 
between d 3y 2_ r 2 and d 3x 2_ r 2, because the overlaps between these two sets of wave functions 
are larger than 97%. This result is in good agreement with the orbital order measured in 



LaMnOs using resonant X-ray scattering [22 . By substituting the orbital wave functions 



into the Hamiltonian Eq. (2), it is straightforward to obtain an anisotropic DE Hamiltonian 
for the charge carriers 



H eff = ~tJ2 W ( X + S ' ' S i)/ 2 d i d 3 ' ( 5 ) 



i.) 



where gr^ assumes different values g» and g± in the ah plane and along the c axis. For 
the antiparallel orbital order, g\\ = 1/2 and g± = 1, the motion of the charge carriers is 
three-dimensional. For the parallel orbital order, g» =3/2 and g± = 0, the charge motion 
is of two-dimensional nature. According to Eq. (5), the charge band width is W = 8t 
for the antiparallel orbital order, and W = 12t for the parallel orbital order, indicating 
that the kinetic energy term favors parallel orbital order. The kinetic energy will compete 
with the effective spin-orbital coupling which is favorable for antiparallel orbital order. The 
spin-orbital coupling is proportional to riiUj ~ (1 — x) 2 according to Eq. (2). For small 
x, the spin-orbital coupling dominates and leads to antiparallel orbital order, in which the 
charges can move in three dimensions. This explains the appearance of the FM state for 
0.1 ~ x ~ 0.4 in the phase diagram Fig. 2. For large x, the spin-orbital coupling reduces, 
the kinetic energy dominates and leads to parallel orbital order. In this case, the charge 
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hopping is confined in two-dimensional planes, and so is the DE ferromagnetic coupling. The 
net interplane spin interaction then comes only from the antiferromagnetic superexchange 
coupling. As a result the A-AF order occurs. 

In summary, we have established a modified DE model to describe the spin and orbital 
properties in the Mn oxides. The obtained phase diagram is consistent with experimen- 
tal observations in the Mn oxides with relatively large conduction band width, such as 
Pri-^Sr-rMnOs and Lai-^Sr^MnOs, where the effect of long-range Coulomb interactions 
may be unimportant and the charge ordering is absent or rather weak. We expect that the 
present model can be further extended to describe the CE-AF state with charge ordering ob- 
served in systems with relatively narrow conduction band by including long-range Coulomb 
interactions. 

This work was supported by the Texas Center for Superconductivity at the University 
of Houston through a grant from the state of Texas, by Texas ARP No: 3652707, and by 
the Robert A. Welch foundation. 
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FIGURES 



FIG. 1. Phase diagram at x = in the A vs J'/ J plane where J = J s (l + 7). 

FIG. 2. Phase diagram in the ksT/t vs x plane with PM as the paramagnetic phase. The 
dashed line indicates the boundaries of the orbital phases. Inset is the orbital phase diagram, 
where "antiparallel" stands for alternating orbital order between (d 3z 2_ r 2 + d x 2_ y 2)j \[2 and 
(d 3z 2_ r 2 — d x 2_ y 2) I 'y/2, and "parallel" indicates uniform d x 2_ y 2 orbital order. 

FIG. 3. Calculated magnetization curves for several values of doping concentration x. 

FIG. 4. Calculated magnetization curves for x = 0.5 at several values of magnetic field. 

FIG. 5. Illustration of the spatial electron density distribution in (a) alternating orbital ordering 
state between (d 3z 2_ r 2 + d x 2_ y 2)j\f 7 l and (d3 2 2_ r 2 — d x 2_ y 2)/\^2, and (b) uniform d x 2_ y 2 orbital 
ordering state. 
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Fig. 1 (Sheng etal.) 




Fig. 2 (Sheng etal.) 
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Fig. 3 (Sheng etal.) 




Fig. 4 (Sheng etal.) 




Fig. 5 (Sheng etal.) 



